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^ Abstract 

^ We investigate the fractal structure of 2d quantum gravity coupled to matter by 

cn 

measuring the distributions of so-called baby universes. We demonstrate that the 
method works well as long as c < 1. For c > 1 it is not clear what distribution 
to expect. However, we observe strikingly similar distributions for various kinds of 
matter fields with the same c. This indicate that there might be some range of c > 1 
where the central charge of the matter fields alone determines the fractal structure 
of gravity coupled to matter. The hypothesis that the string susceptibility 7 = 1/3 
is found to be compatible with the data for 1 < c < 4. 
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1 Introduction 



The fractal and self-similar structure of 2d quantum gravity depends on the central 
charge of the matter fields coupled to gravity. For c < 1 the relation is given by the 
famous KPZ formula [1] 



c - 1 - J 1 - c 25 - c 
A ; 

' 12 ' ^ ^ 

where 7 denotes the entropy (or string susceptibility) exponent of the theory. The 
detailed relationship of 7 to the fractal structure of the ensemble of random sur- 
faces dictated by the partition function was analyzed in a recent paper by Jain and 
Mathur. 7 determines the distribution of so-called baby universes which can branch 
off from the "parent universe". In [3] it was shown that the measurement of baby 
universe distributions is an efficient way to determine 7 from numerical simulations. 
In this letter we report the results from extensive simulations using this method to 
determine the fractal structure of 2d gravity coupled to matter fields. 



2 The method 

Let us briefiy summarize the method: We discretize the theory of 2d gravity coupled 
to matter by the use of dynamical triangulations [4, 5, 6]. The partition function 
can be written as 

Z(/.,/3) = ^Z(iV,/3)e-^^, (2) 

N 

where N denotes the number of triangles, fi the cosmological constant and /3 some 
matter field coupling. Critical behaviour for // — )• /ic is obtained if Z[N^ fi) has the 
following asymptotic behaviour for large N: 

Z{N, fi) ~ e^^^f^^^ N^m-^ \og»(P) iv. (3) 

We have in (3) allowed for a logarithmic correction term of the type that is known 
to play a role for c = 1. 

We assume spherical topology and define for a given triangulation a (minimum 
neck) baby universe as follows: If we cut open the triangulated surface consisting of 
N triangles along a closed loop made of three different links it will be separated in 
two disconnected parts, consisting of B and N — B triangles, where B < N/2. We 
call the part of triangulation containing the B triangles the baby universe. From 
(3) it follows that the average number of baby universes consisting of B triangles on 



2 



a closed triangulated surface of spherical topology made of N triangles is given by 

hN{B) ~ [{N - 3)3]^^^^-^ [log(iV - B) log BY^^\ (4) 

It is now clear that a measurement of n]\j{B) for a fixed large N in principle should 
be sufficient to determine both 7(/3) and 

The matter fields we couple to 2d gravity will be a variety of (multiple) q-siaie 
Potts models, q = 2,3,4, all with < c < 4. The matter fields are placed in the 
center of the triangles and the implementation is standard (for details we refer for 
instance to [7]). In the case of the q-siaie models we have a coupling constant fi 
multiplying the matter part of the action. If we have a multiple q-siaie model we 
still consider only the case where a single common coupling constant multiplies the 
different identical copies in the action. We used canonical Monte Carlo simulations, 
and lattice sizes ranging from 1000 to 2000 triangles. To update the spin configu- 
rations a Swendson-Wang cluster algorithm was used, while for the graphs we used 
a standard link-fiip algorithm. Each Monte Carlo sweep consisted of one cluster 
update and N link fiips. Usually 2.5 X 10^ sweeps where performed, but close to fic 
we used up to 3 X 10'' sweeps. 

For c < I we know that the models in the infinite volume limit (// = /ic) un- 
dergo a phase transition from a phase with no magnetization to one with constant 
magnetization for a certain value fic of fi. At this point 7(/3) changes from its value 
without coupling to gravity to a new value, given by (I). This gives us non-trivial 
tests of our method. On a finite lattice the jump from the pure gravity cannot take 
place at a single fi value, and we see a peak around the critical value. However, as 
will be clear from the numerical results, the peak is very sharp. 

For c > I we have no analytical results. Numerical studies suggest that we 
will still observe a transition between a phase with zero total magnetization and a 
magnetized phase. However, as reported in [7] the picture seems somewhat more 
complicated if c is large since the back-reaction of the matter fields on gravity was 
largest before the transition, i.e. for fi < fic- We will return to a discussion of these 
results in sec. 4. 

3 The measurements 

We have measured njv(-B) with high statistics and for a wide range of couplings in 
the following models: single q = 2,3,4-state Potts model (c = 1/2, 4/5 and I), a 
double q = 3-state Potts model (c = 8/5) and a quadruple q = 3-state Potts model 
(c = 16/5). We also looked at double q = 4-state and quadruple q = 2-state models 
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(c = 2). and a quadruple q = 4-state model (c = 4), for couplings around fic- The 
measured distributions n]\j{B) are fitted to 



lognAr(5) = A + (7-2)log 



Bil - ^] 



+ alog[logBlog{N -B)] + ^. (5) 



N'l " B 

The last term is a "phenomenological" finite size correction term. The asymptotic 
form (3) of Z{N,fi) is only assumed valid for B » 1. A simple first correction 
could be of the form 

It leads to the last term in (5), and has been shown to improve the results for small 
B considerable [3]. 

In fig. I we show the results of the fits for 7, with and without the additional 
parameter a from the logarithmic corrections, for values of fi near the critical value 
fic for the various models. The inclusion of a in the fit increases the errorbars on 7 
for the obvious reason that it can be hard to distinguish the two contributions over 
a limited interval of B. 

The general picture is as expected: we see a peak in 7(/3) near /3c, while 7(/3) 
has essentially reached the c = value 7(c = 0) = —1/2 as soon as |/3 — fic\ > 0.1. 
This is a clear indication that 7 changes at the critical point. Naively one would 
expect the peak to get more narrow if we increase the number of triangles N used in 
the simulation, but we have only seen a weak dependence on N. If we compare the 
functions 7(/3) for c < 1 with the corresponding function for c = 1 we see a marked 
difference in the shape. For c < 1 there is a very narrow and pronounced top in the 
peak as already mentioned, while the peak is rather broad for c = 1. Further the peak 
value of 7(/3) gives the correct result for c < 1 when we do not include logarithmic 
corrections. If we include logarithmic corrections the value will be somewhat to high. 
Since we know from analytical results that there are no logarithmic corrections for 
c < 1 we conclude that it is difficult to include a as a free parameter in the fit. 
For c = 1 the situation is different. The fit without logarithmic correction gives 
too small a peak value for 7(/3). If fact the value is smaller than the corresponding 
values for c = 1/2 and c = 4/5. The very sharp peaks of the c < 1 curves are simply 
missing. Including the logarithmic correction on the other hand gives the expected 
value of 7 at the peak (7 = 0.05 it 0.08 for N = 1000). The logarithmic exponent 
was harder to determine but we got a = —1.6 it 0.3 (compared to the expected value 
of a{/3,) = -2). 

Our conclusion is that the baby universe technique for extracting 7 works very 
well for c < 1, and that triangulations of the size 1000-2000 triangles and baby 
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Figure 1: Measured value of 7. In all figures circles denote lattice size N = f 000 and 
squares N = 2000. In fig. fa we shown f Ising model (open symbols) and 1 q = ?>- 
state Potts model (filled symbols). The fits are without logarithmic corrections. In 
fig. fb we show 1 q = 4-state Potts model and in fig. Ic2 q = 3-state Potts models. 
Open symbols show fits without logarithmic correction and filled symbols fits with 
logarithmic corrections. 
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0.64(9) 


-3.16(92) 


180 
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5 Gaussian 


0.217(10) 


1374 


0.97(17) 


-4.34(75) 


155 
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Table 1: Fxponents and values for fits to lognjv(-B) for c > le with i? > 11 and 
245 data points. The following functional forms are used: 

(a) A + (7 - 2) log [B{1 - BjN)] + CjB 

(b) A + (7 - 2) log [B{1 - BIN)] + a log [log51og(iV - B)] + C/B 

(c) A + (1/3 - 2) log [B{1 - BIN)] + a log [log51og(iV - B)] + CjB 



universe distributions of a size from 10-500 are sufficient to reflect in a precise way 
the fractal structure of 2d quantum gravity coupled to matter with c < 1,. However, 
this conclusion is based on the knowledge that there are no logarithmic correction 
for c < 1. 

Let us now turn to a discussion of the results for c > 1. It has often been argued 
that c > 1 has a lot in common with gravity in higher dimensions. The program of 
dynamical triangulation has been formulated and implemented in higher dimensions 
[9], and one of the lessons is that we cannot take the functional form (3) for granted. 
It seems that the generic form is rather like [10] 

Z(iV) ~ e^=^-^^^ a<l. (7) 

We have tested whether Z{N, fi) for c > 1 should fit better to such a functional 
form than to (3), but found that it is not the case. Fven though in some cases 
the of the fits where not bad, the fitted parameters had absurd values. As a 
consequence we will assume that (3) is a reasonable functional form. However, for 
c > 1 we have no a priori knowledge of the logarithmic corrections. It is notoriously 
difficult to distinguish between logarithmic and small power law corrections, but 
the good results for the c = 1 case might give some hope that one can extract a 
reliable the value for 7. As will be obvious from table 1 and sec. 4 this will not 
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be the case^. The trend of broader peaks and larger values of a if extracted from 
(3) continues with increasing c. Again one could be worried that we use too small 
triangulations. To test this we have performed a simulation with two q = 'i models 
coupled to gravity and N = 100.000. This allowed us to measure baby universes 
with reasonable statistics up to a size of 5000, i.e. a factor 10 more for the smaller 
surfaces (N= 1000 and N= 2000) used above. We saw however no drastic change in 
the extracted value of 7. If extracted with a set to zero, 7 changed from —0.212(15) 
to —0.175(24). If extracted with logarithmic correction the numbers were 0.11(7) 
and 0.19(9), respectively. In both cases the of the fits were acceptable and we 
have no numerical reason for preferring one kind of ansatz from the other. 

The results obtained for c > 1 are shown in table 1 together with results ob- 
tained for Gaussian models [3]. For the same c we get the same 7's and, within the 
precision available, the same a's^. One is lead to the hypothesis that there might 
be universality even for c > 1, at least in some region 1 < c < Cq (the largest value 
of c where we have compared distributions in detail is c = 4). 

This becomes even more striking when one makes a direct comparison between 
the distributions njv(-B) of baby universes for c = 1,2 and 4. In fig. 2 we show 
the measured distributions (normalized by the the distribution for pure gravity) 
for the following models: For c = 1 1 Gaussian and 1 q = 4-state Potts model, 
for c = 2 2 Gaussian, 2 q = 4-state Potts and 4 Ising models, and for c = 4 4 
Gaussian and A q = 4-state Potts models. The advantage of this representation 
is that it is independent of any fits and assumptions about the functional form of 
the distributions. We see a pronounced universality in the curves for fixed central 
charge. It seems as if the back-reaction of matter on gravity only depends on the 
central charge even if c > 1. 

In order to make the claim more substantial we have statistically tested the hy- 
pothesis that the distributions are identical, module a possible non-universal normal- 
ization and phenomenological correction terms. So we fitted the difference between 
two distributions and ^2^^, with the same c, to the form 

logni^(5) - logn2^(5) = oi + ^, (8) 

and tested the goodness of the fits. The fits where made with i? > 11 and 245 
data points. In table 2 we show the resulting for the fits and the corresponding 

similar conclusion as been reached in two recent papers using different observables [11, 12]. 
^By 7 and a for the spin models we mean the values measured at the (finite size) critical /3 as 
identified by the peak value o{j(p). The Gaussian models are alway critical (before coupling to 2d 
gravity) and any coupling constant can be absorbed in a redefinition of the cosmological constant. 
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Figure 2: Normalized distributions of baby universes for various matter fields cou- 
pled to 2d gravity. 
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Table 2: Quality ol the fits to (8) lor comparing distributions ol baby universes with 
same central charge 

significance levels -P(x243 — X^)- can be seen the significance levels are in all cases 
larger than 50%, so we conclude that the data is consistent with the hypotheses. 

Contrary to the observables used in [7, 8], where Gaussian models and multi- 
ple Potts models coupled to 2d gravity were also compared, and which showed a 
surprising universality between various Potts models with the same central charge, 
but some deviation between the Potts models and Gaussian models, the distribu- 
tions of baby universes should (for large B) be characterized only by their critical 
exponents. They are consequently much better observables when it comes to a 
comparison between different models. 

4 The quest for a c > 1 hypothesis 

The results reported above show that distributions of baby universes provide us 
with a convenient tool for analyzing the interaction between matter and 2d gravity. 
They also suggest that the fractal structure of the random surfaces depends only 
on the central charge of matter, at least in the range 1 < c < 4 we have tested 
by numerical simulations. Unfortunately the method is unable, with the present 
quality of the numerical results, to nail down the correct values of 7 and a in (3) 
as they correspond to terms with too similar functional behavior and a change in 
one can be absorbed in the other without seriously affecting the quality of the fit. 
From this point of view it would be very interesting to have an educated guess of 
a more restricted functional form than the one given by (3). Let us end by a brief 
discussion of some recent speculations in that direction [13, 14, 15]. 

It has been argued that 7 > implies 7 = l/{n -|- 1), n > 1, [13, 14]. Specific 
matrix models which realize 7 = l/(n-|-l) are known [17]. For spin systems it 
is natural to imagine that 7 = 1/3 is realized, at least for multiple spin systems 
with c >> 1. The argument goes as follows: The value 7 = l/{n -\- 1) refiects a 
fractal structure of the surfaces where the surface is glued together of "bubbles" of 
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Figure 3: The logarithmic correction exponent a lor the fits (5) and (10). 



7 = —1/n matter. The results in [7] strongly suggest that the individual bubbles 
are already magnetized in some range of /3 < /3c and that the phase transition to a 
magnetized phase is just an alignment of the spins in the individual bubbles. This 
would imply 7 = — 1/2 (pure gravity) and therefore 

7 = 1- (9) 

An identical scenario was obtained independently by Wexler using mean field theory 
[15]. While this result might well be correct for multiple Potts models with c >> 1 
it is not clear that it is valid for Gaussian models for c >> 1, where other mean 
field results suggest that 7 = 1/2 or 7 = — c/2, depending on the weight used for 
the triangulations [16]. However, it is tempting to conjecture that 7 = 1/3 should 
be valid all the way down to c = 1 for multiple Potts models, and since we have seen 
that the baby universe distributions (within the numerical accuracy) for small c > 1 
are the same for Gaussian and Potts models, there should be a range 1 < c < Cq 
where the Gaussian models also have 7 = 1/3. If we combine this conjecture with 
the numerical observation that the baby universe distributions change continuously 
in some range c > 1, we are lead to the suggestion 

Z,{N) = log"(') iV, c> 1, a{c) -00 for c ^ 1. (10) 

We have tested whether this conjecture is in accordance with our baby universe 
measurements for c > 1. The result is shown in table 1, where we show extracted 
values of a and the quality of the fits. We see that of fits made with the assump- 
tion 7 = 1/3 and with 7 as a free parameter is comparable. This illustrates further 
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that additional information is needed if we want to have some hope of extracting 
reliable values for 7. It also shows that the conjecture (10) cannot be ruled out. 
It is interesting to note that we get a very different functional dependence a[c) for 
the logarithmic exponent a when we compare the results where 7 is fixed to 1/3 
and results where 7 is arbitrary. This is shown in fig. 3, and we observe that the 
functional form a[c) for 7 = 1/3 indeed is compatible with (10), except that we 
cannot expect a[c) — )• —00 for finite size lattices. 
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